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Theory of Sodium Ordering in Na x Co02 
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The ordering of Na ions in Na^OoC^ is investigated systematically by combining detailed density 
functional theory (DFT) studies with model calculations. Various ground state ordering patterns 
are identified, and they are in excellent agreement with available experimental results. Our results 
suggest that the primary driving force for the Na ordering is the screened Coulomb interaction 
among Na ions. Possible effects of the Na ordering on the electronic structure of the C0O2 layer are 
discussed. We propose that the nonexistence of a charge ordered insulating state at x = 2/3 is due to 
the lack of a commensurate Na ordering pattern, whereas an extremely stable Na ordering at x = 0.5 
enhances the charge ordering tendency, resulting in an insulating state as observed experimentally. 

PACS numbers: 61.50.Ah,61.66.-f,61.18.-j 
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The recently renewed research interest in Naj;Co02 
since the discovery of superconductivity in the hy- 
drated materials has revealed a range of interesting and 
intriguing properties of this system. One of the most 
interesting discoveries is the determination of its phase 
diagram as the doping level x is varied Q: Two metal- 
lic phases at low and high dopings are separated by an 
insulating state at x = 0.5. Along with this insulating 
state there is strong evidence of an ordered Na layer and 
commensurate charge orderings in the C0O2 layer. This 
raises the possibility of a subtle interplay between the Na 
ordering and the charge ordering in the Co0 2 layer[ll3. 
At certain Na compositions, this interplay might have 
profound effects on the electronic or magnetic properties 
of the C0O2 layer if a particularly stable Na ordering pat- 
tern exists and whether a commensurate charge ordering 
is allowed. Therefore, a unified theory of the Na order- 
ing mechanism may be critical for a better understanding 
of the electronic structure in this system. The observed 
Na ordering pattern ^aOi at x = 0.5 agrees well with our 
theoretical prediction [5j , and this motivates us to extend 
our investigation on Na ordering at other compositions 
studied by recent electron diffraction experiments [j] . In 
this paper, we identify ground state Na ordering patterns 
at various compositions and show that the Na ordering is 
primarily driven by the intra-plane screened electrostatic 
interactions. The details of the electronic structure of the 
C0O2 plane, on the other hand, have minimal effects. 

The compound Na^CoO^ assumes a layered structure 
which consists of alternating triangular C0O2 and Na 
planes. Although there have been reports d, S S @ of 
several stacking patterns, we shall restrict our discussion 
to the so-called 7 phase, which is usually observed for 
low to intermediate Na compositions. There are two dis- 
tinct Na sites within a given plane, denoted as Na(l) and 
Na(2) (Wyckoff indices 2b and 2d). The Na(l) site, being 
directly between two Co ions, is slightly higher in energy 
than the Na(2) site@. Therefore, Na ions will normally 
prefer occupying the Na(2) sites. 



We calculate the total energy of 82 ordered structures 
of Naa;Co02 using first principles techniques. The choice 
of structures for the ground state search was guided 
by the MAPS code [T(|- Our calculations are based 
on density-functional theory (DFT) 0, within the 
local density approximation (LDA) with the Ceperly- 
Alder exchange-correlation functional [ l3l Il4 | . Ab initio 
Troullier-Martins pseudopotentials [15| are used. Calcu- 
lations are performed using the SIESTA code [3 G3, 
which expands the Kohn-Sham wavefunctions in a linear 
combination of atomic orbitals (LCAO). A double-zeta 
plus polarization (DZP) basis is used. The energy cutoff 
for the charge density grid is 1000 Ry. The irreducible 
Brillouin Zone is sampled using Monkhorst-Pack grids 
|18| with a density of k-points equivalent to 12 x 12 x 2 for 
a primitive unit cell. Atomic positions are relaxed until 
forces are smaller than 0.005 eV/A. We keep the lattice 
constants fixed at a = 2.82A and c = 10.89A. Although 
the lattice constants of Na.j;Co02 will vary slightly with 
Na concentration, our results are not sensitive to this 
small variation. 

The energetics of Na ordering in Naa;Co02 can be stud- 
ied using the well-stablished formalism for binary alloys 
|19|. The fundamental quantity that determines the sta- 
bility of a given arrangement r\ of Na ions and vacancies 
with Na composition x is the formation energy AE 11 : 



AEl FT = E r > - x E NaCo ° 2 



~ (l-x)E 



C0O2 



(1) 



where E v , E NaCo ° 2 and E Co ° 2 are the total energies 
per site of the structure 77, and of the pure compounds 
NaCo02 and C0O2, respectively, calculated within DFT. 
Fig. n shows the calculated formation energies of all 82 
structures (open dots). We highlight the ground state 
structures (red dots) joined together by the hull-shaped 
curve (red line). Structures with formation energy above 
the convex hull are unstable against phase separation into 
the two ground states with nearby compositions. Red ar- 
rows indicate the structures proposed by Zandbergen et 
al. (at x = 0.11, x = 0.30, x = 0.33, x = 0.50, 
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FIG. 1: DFT formation energies for 82 ordered structures 
of Nai,Co02. The red dots represent the ground-state struc- 
tures, joined together by the convex hull (red line). The ar- 
rows indicate the experimentally observed ground states Q- 



x = 0.71 and x — 0.75) based on electron diffraction 
experiments. The agreement between theory and ex- 
periment is remarkable in that all experimentally pro- 
posed structures are indeed ground states within DFT. 
Although the structure at x = 0.30 appears to be an ex- 
ception (it is not a ground state within DFT), it is too 
close to the convex hull to be neglected as such, given 
the uncertainties in the calculation. Moreover, we find 
additional ground states at x — 0.20, x = 0.25, x = 0.60, 
x = 0.80 and x = 0.83. A structure at x = 0.40 is also too 
close to the convex hull to be ruled out. These structures 
have not been observed experimentally, perhaps because 
samples at these particular compositions were not an- 
alyzed. An ordered structure at x = 0.64 also appears 
experimentally 4]. According to our results, at this com- 
position the sample should be phase-separated between 
the x — 0.60 and the x = 0.71 ground states. All ground 
state structures and the near ground states at x = 0.30 
and x = 0.40 are depicted in Fig|21 We also display the 
lowest energy structure at x = 0.67 (not a ground state), 
for discussion purposes. 

Before discussing in detail the various ordered struc- 
tures, let us first address the driving force for Na order- 
ing in this system. The simplest explanation would be 
electrostatics since the Na atoms lose one electron each 
to the C0O2 network, they arrange themselves so as to 
minimize electrostatic energy. We perform model calcu- 
lations to test this hypothesis, in which the Ewald (ionic) 
energies of the different arrangements of Na +1 ions are 
computed. In these calculations, we assume that the elec- 
tronic charge is spread over the C0O2 planes neighbor- 
ing the Na-vacancy layer, which is consistent with our 
previous calculations 0. Also, the chemical difference 
between Na(l) and Na(2) sites is simulated by adding an 
energy t\ per Na ion to the Ewald energies for a Na(l) 
site. Therefore, the total energy per unit cell for this 



(a),v = (I))* = 0.11 (c) x = 0.20 




FIG. 2: Relevant ordered structures of Na^CoO 2. The dashed 
lines represent the projected triangular lattice of Co atoms. 
Large red open circles and small red dots represent projected 
Na positions in different planes, at z = and z — 0.5c, re- 
spectively. The Na(l) sites fall on top of the triangular lattice 
sites and Na(2) sites fall in the center of the triangles. Blue 
lines indicate the minimum unit cell in each case. A detailed 
discussion of each structure is presented in the text. 

electrostatic model is defined as: 

El odel ^El wald + N^e 1 , (2) 

where N± is the number of Na atoms in Na(l) sites for 
the structure rj. 

Figure |21 shows, for each structure, /S.E V DFT (Eq.Q) 
versus the calculated formation energy using the above 
model (A£^ ode/ ). We clearly see that these two quan- 
tities are proportional: A£^, FT = nAE^^. The best 
fit (rms deviation of 0.016 eV) is obtained for k = 0.514 
and an on-site energy difference between Na(l) and Na(2) 
sites Kei = 67 meV. A value of k smaller than 1 is ex- 
pected because of screening and partial ionization effects, 
which are absent from the model calculations and natu- 
rally included in the DFT calculations. Therefore, DFT 
formation energies can be described with good accuracy 
by a model of screened electrostatic interactions, show- 
ing unambiguously that the dominant driving force for 
Na ordering in Na x C0O2 is electrostatic. 

We now discuss the geometries of the structures shown 
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FIG. 3: DFT versus electrostatic model formation energies. 

in Fig. |21 In the low concentration regime, we find a 
ground state at x = 0.11 (1/9) (Fig. Hfb)), correspond- 
ing to an in-plane arrangement of Na "impurities" in a 
3x3 triangular lattice. Surprisingly, the lowest energy 
structure corresponds to placing both Na atoms (one in 
each plane) in the less favorable Na(l) sites. This struc- 
ture prevails with respect to other arrangements because 
it minimizes the inter-plane electrostatic energy: The Na 
atoms can be as far as possible from each other. Inter- 
estingly, this arrangement corresponds to placing each 
Na atom in the middle of the projections of the three Na 
atoms of neighboring planes, and it is consistent with the 
experimental observation of a y/3 x V3 superstructure in 
the diffraction patterns of Ref. Q if a random stacking 
of such arrangements occurs. 

At x = 0.20, the ground state has a v7 x V7 super- 
cell, shown in Fig.|2Jc), with both Na ions in Na(2) sites. 
Periodic repetition of this cell leads to an arrangement 
of Na ions along rows perpendicular to the [110] direc- 
tion (the horizontal direction in Fig. 2). The distance 
between Na ions in a row is y/3a. Such pattern of rows 
have been experimentally observed for different composi- 
tions and have been proposed to be the guiding principle 
for Na ordering :4] . Indeed, we observe this pattern in all 
other Na compositions except x = 1/9, 1/4, and 2/3 (not 
a ground state) . The same pattern of rows is observed in 
the lowest energy structures calculated using the simple 
electrostatic model. Therefore, the rows arise naturally 
from the minimization of electrostatic energy, combined 
with the constraint that only discrete lattice sites can be 
occupied by Na ions. 

Fig. Hfd) shows the ground state for x = 0.25, a 2 x 2 
triangular lattice. In this case, both Na ions occupy the 
more preferable Na(2) sites. In Fig. |2fe) , the lowest en- 
ergy structure for x = 0.30 is shown. It has precisely 
the same in-plane arrangement proposed by Zandbergen 
et aZ.Q. Again, it corresponds to rows of Na ions per- 
pendicular to [110] and all ions occupy Na(2) sites. As 
discussed above, it is not a ground state in our calcula- 
tions (by merely 1 meV/site). 



At x = 0.33 (1/3), the arrangement is again a trian- 
gular lattice (V3 x y/3), as shown in Fig. Elf)- This 
structure is exactly the same as proposed by Zandbergen 
et al. 4]. The x = 1/3 concentration is of particu- 
lar interest due to the observed superconductivity in the 
hydrated system at x ~ 0.3. Although the role water 
plays is still unknown, the screening of the Na potential 
by water molecules might suppress the tendency of de- 
veloping a commensurate charge ordering pattern in the 
C0O2 layer, which may then lead to a more homogeneous 
electronic system and favor a superconducting state over 
other competing phases at low temperature. Possible 
charge disproportionation and gap opening at this dop- 
ing level in the unhydrated system has been studied pre- 
viously |2(j . However, so far such an insulating state has 
not been observed. 

Fig. EJg) shows the lowest energy arrangement for x — 
0.40 (not a ground state, by 0.6 meV/site). It also follows 
the rows pattern, and the ratio between Na(l) and Na(2) 
occupancies is 1. 

In Fig. |2fh), the ground state for x = 0.50 is shown, 
which is a v3 x 2 supercell. The arrangement also fol- 
lows the pattern of rows and it has also been predicted 
theoretically |5j and observed experimentally [j]. This 
structure is particularly stable since all other in-plane 
arrangements with a similar unit cell at this composi- 
tion are at least 25 meV/site (or 200 meV/cell) higher 
in energy. Note again that the ratio between Na(l) and 
Na(2) occupancies is 1. This ratio tends to be higher for 
structures with intermediate composition, since the on- 
site chemical energy difference between Na(l) and Na(2) 
sites becomes less important than arranging the Na ions 
as far as possible to minimize electrostatic energy. The 
existence of a particularly stable Na ordering pattern, to- 
gether with the observed insulating state at this compo- 
sition [2j, strongly indicate an interplay between the Na 
ordering and the charge ordering in the C0O2 layer in this 
system. Our results suggest that, if such an interplay ex- 
ists, it is more likely that an ordered Na potential drives 
the charge ordering or enhances this tendency. Studies 
of charge ordering in this system with the presence of a 
commensurate Na ordering will be reported in a separate 
paper. 

At x — 0.60, the ground state structure is again a 
V7 x v7 supercell, shown in Fig. |2fi) with a pattern 
of rows and a Na(l)/Na(2) occupancy ratio of 1/2. The 
composition x = 0.67 (2/3), like x = 1/3, has been the 
subject of many discussions due to its ap parent p roximity 
to charge ordering in a triangular lattice [20L l2ll |22| . The 
lowest energy structure corresponds to 3 x 1 supercell, 
shown in Fig. |2[j). This structure is clearly lower in 
energy than any honeycomb arrangement of Na ions. It 
is obvious that a honeycomb arrangement would be more 
compatible with the proposed charge ordering natternpol 
0,|22]. Therefore, we propose that the nonexistence of a 
charge-ordered insulating state at this composition is due 
to the lack of a stable and commensurate Na ordering. 

The ground state structure at x = 0.71 (5/7) is shown 
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in Fig. Oik). It is impressive that such a complex ar- 
rangement with a large supercell (\/l3 x \/l3) corre- 
sponds to a ground state structure, as seen experimen- 
tally p|. For this structure we have a Na(l)/Na(2) oc- 
cupancy ratio of 2/3 and once more the rows pattern is 
followed. In the high x limit, this pattern becomes sim- 
pler to analyze: Low energy structures are composed of 
rows of Na(2) ions intercalated by rows of vacancies and 
Na(l) ions. For instance, for the x — 0.71 structure the 
repetition is 3xNa(2)-vac-2xNa(l)-vac. 

For x = 0.75 we find again the same ground state as 
proposed by experiments .4] , a V3 X 4 supercell shown in 
Fig. in this case, the pattern of rows is 3xNa(2)- 

vac-3xNa(l)-vac, i.e., there are equal numbers of Na(2) 
and Na(l) ions. This structure is lower in energy than 
that proposed by Shi et fl/.|23|. Measurements of local 
distortions of Na(2)C>6 polyhedra at x = 0.75 have been 
recently performed by neutron diffraction |24| . Appar- 
ently, Na(2) ions move to off-center positions, resulting 
in two long (2.56 A) and four short (2.32 A) Na-0 bonds. 
Our calculated Na-0 bond lengths for this structure are 
2.57 A and 2.29 A, in excellent agreement with experi- 
ment. 

At x = 0.80 and x = 0.83 (5/6) (Figs. (m) and 
(n), respectively) we find ground-state structures that 
belong to the same "family" as the x = 0.75 structure, 
i.e., they follow similar row patterns with nxNa(2)-vac- 
3xNa(l)-vac rows in a x (n + 5)/2 cell, with n = 3 
for x = 0.75, n = 5 for x = 0.80 and n = 7 for x = 0.83. 
Therefore, it is very likely that beyond x — 0.75 one can 
construct infinitely many ground states (polytipoids) by 
simply intercalating odd numbers of Na(2) rows with the 
sequence of rows -vac-3xNa(l)-vac-. In fact, if one takes 
into account entropic considerations, it is plausible that 
the thermodynamically stable configuration for x > 0.75 
will consist of disordered arrangements of such rows. This 
result may be related to the order-disorder transition re- 
ported at x = 0.75 from neutron diffraction experiments 
Uliil]. Finally for x = 1.0 all Na(2) sites are occupied 



(Fig. Ho)). 

In conclusion, we have carried out detailed DFT 
and model studies on the Na ordering mechanism in 
Na^CoC^. The ordering pattern is non-trivial and is 
sensitive to the Na concentration. The theoretically de- 
termined ordering patterns at various Na compositions 
agree well with those observed in recent electron diffrac- 
tion experiments 4] . In addition, we identify several 
ground states that have not yet been observed. Our re- 
sults indicate that the primary driving force for the Na 
ordering is the screened electrostatic interactions. De- 
tailed electronic structure on the C0O2 layer, on the other 
hand, plays no important role on the ordering of Na ions 
since the coupling between the C0O2 plane and the Na 
plane is much weaker than the intra-plane Na-Na interac- 
tions. The converse effects, however, could be significant 
since the charge and magnetic orderings on the C0O2 
plane happen at extremely low energy scale Q and may 
be vulnerable to external perturbations. Therefore, it is 
likely that the particularly stable Na ordering structure 
at x = 0.5 enhances the charge ordering tendency in the 
C0O2 layer, resulting in a charge ordered insulating state 
at low temperature as observed experimentally 2]. Con- 
versely, the nonexistence of such an insulating state at 
x = 2/3 may be due to the absence of an ordered honey- 
comb arrangement of Na ions at this composition. 
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